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Copyright © 1975, Texas Instruments Incorporated 


The program material contained herein is supplied without 
representation or warranty of any kind. Texas Instruments therefore 
assumes no responsibility and shall have no liability, consequential or 
otherwise, of any kind arising from use of this program material or any 
part thereof. 


ADDENDUM 


Dear Customer: 
Please mark the following corrections in the Statistics Library Program 
Manual ST1: 


PROGRAM PAGE SHOWN AS CHANGE TO 
ST 1-03-1 15 196 47 *1' 196 87 *1’ 
ST 1-05 25 218 18 C’ 216 36)" 
ST 1-06 29 120 04 4 129 04 4 
ST 1-15 re 132 04 4 133 04 4 
ST 1-25 103 002 02 INV 002 22 INV 
ST 1-14 66-67 (see other side) 


DO NOT RERECORD THE MAGNETIC CARDS. The prerecorded magnetic 
cards have been correctly recorded at the factory. 
Change the first equation on page 18 to read: 
U;,, = 79(U; = 10'°)(Mod 10°) 
On page 65, the minus sign should be deleted from the display column 
in Step 6. 
‘Note 2 on page 68 should be changed to read |X| < 21.097 
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ADDENDUM 
SR-52 Program Manual ST1 Statistics Library 


This manual was printed with an incorrect program listing for program 
S1T1-14 Histogram Construction. DO NOT rerecord the magnetic card for 
this program. The correct program listing follows. 


Sibel PROGRAM LISTING 


war [ver [own [ev | cena [ev | oan | ver 


01 

0S. 5 

95 = 

80 *if pos 





*Denotes 2nd function key 


ST1-14 


PROGRAM LISTING 


86 “rset 


*LBL 





“Denotes 2nd function key 
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ST1-01. BASIC STATISTICS FOR ONE OR TWO VARIABLES. ... . 
Calculates the means, standard deviations, and standard 
errors for one or two variables using both the normalized 

~ and unnormalized methods. Also calculates the covariance 


' and correlation coefficient for the two variables. 


ST1-02. PERMUTATIONS AND COMBINATIONS .......... 
Calculates the permutation and combination for 
agivennandr. 


57 1-03-1,2. MEANS AND MOMENT) 6. bi bw me CONT Hi oe Pa epi 
The arithmetic, geometric, harmonic, and generalized 
means, the first four moments, and the kurtosis and 
skewness of distribution are calculated for grouped or 
ungrouped data. 


ST1-04. RANDOM NUMBERGENERATOR............. 
Uniformly distributed and normally distributed 
random numbers are generated. | 


ST1-05. ONE WAY ANALYSIS OF VARIANCE. ........... 
Performs a one way analysis of variance on k treatment 
groups. Includes the F statistic, degrees of freedom, 
and sums of squares. 


ST1-06. t STATISTIC EVALUATION (PAIRED OBSERVATIONS) 
. Uses t statistics to test the difference between the means 
-of two normally distributed populations. 


ST1-07.t STATISTIC EVALUATION (TWO SAMPLE TEST). ..... 
Evaluates the t statistic with n, + n, — 2 degrees 
of freedom. 


Be eR PRP CR TRU IN  k s w k C w le a 
The slope, intercept, estimated values for x or y, and 
coefficient of determination for the linear least-squares 
fit of given points (x,y) are calculated. 
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7108. POWER CURVERTT... . . wo Oe es Ss ae 


The slope, intercept, estimated values for x ory, and 
correlation coefficient for a power least-squares 
fit of given points (x,y) are calculated. 


or 1-10, EXPONENTIAG CAIIVE FEE ee a? Sr eee at 


The slope, intercept, estimated values for x or y, and 
correlation coefficient for an exponential least-squares 
fit of given points (x,y) are calculated. 


ST1-11. LOGARITHMIC CURVE FIT .......... eee ae 


The slope, intercept, estimated values for x or y, and 
correlation coefficient for a logarithmic least-squares 
fit of given points (x,y) are calculated. 


ST1-12-1,2. MULTIPLE LINEAR REGRESSION .......... 


The coefficients of the linear least-squares fit of 
given points (x,y,z) are calculated as well as the z 
corresponding to given x and y for this line. 


>! 1-13-1,2. QUADRATIC REGRESSION .. 03k bok a ee ee 


The coefficients of the quadratic least-squares fit of 
given points (x,y) are calculated and the value of y 
for a given x is determined for this curve. 


ST1-14. HISTOGRAM CONSTRUCTION ............. 


Constructs a histogram over a given range and given 
observed data points. The frequencies are calculated for 
each of 12 cells and the mean and variance are calculated 
for the entire range. 


oTT-15. NORMAL DISTRIBUTION. ... ws we ee ce kk ee ? 


The standard normal distribution for a given x is 
calculated using a probability function. 


ST1-16. CHI-SQUARE DISTRIBUTION. .........2.... 


The chi-square density function and probability I'(4) 
are calculated, given the degrees of freedom and x. 


Bee Pe ee I a ie aia hha an cg BL ke VAL se hee wk 


The integral of the t distribution is calculated, 
given the degrees of freedom and x. 
\ 
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818. FiOS Vee Ne WIGS SINE OR SAI Dea ie ai a 
The probability of an event X > x (area of the ‘“‘upper tail’) 
for the F distribution is calculated for two given degrees 
of freedom and x. 


~ $71-19. BIVARIATE NORMAL DISTRIBUTION ........... 84 


The bivariate density function is calculated given the 
means of x and y, standard deviations of x and y, and 
the correlation coefficient between x andy. 


$T1-20. LOGARITHMIC NORMAL DISTRIBUTION ......... 88 
The median, mode, mean, and variance are calculated 
for the logarithmic normal distribution given m and o”. 
Given an x, f(x) is then calculated for this distribution. 


Oe 1. Wer Liter tere FURIE 6 i ee Oe OE ae 90 
The probability density and cumulative probability 
functions are calculated for the Weibull distribution of 
given m, n, and x. Given Q(x), the program also calculates x. 


Bt 1-22, POISSON DISTRIBUTION. 6 cc cleo Pr AW eae +e ers 92 
Calculates the density function and cumulative 
distribution of the Poisson distribution given m and n. 


ST1-20. CINOMIAL Dre TON Se ee eh eee a ee 94 
The binomial density function, the cumulative distribution, 
mean, variance, and probability of at least k successes are 
calculated given n, p, and kK. 


$T1-24. NEGATIVE BINOMIAL DISTRIBUTION. .......... 98 
The mean, variance, negative binomial density function, 
cumulative distribution, and probability of k or less 
successes are calculated given r, p, and k. 


St 1-20. GEOMETAIG DISTRIBU Trini eee sella 6 ea ek a 102 
The variance, probability, cumulative distribution, and 
geometric density function are calculated given u and x. 


S$T1-26. HYPERGEOMETRIC DISTRIBUTION. .......... 104 
The hypergeometric density function, cumulative 
distribution, mean, and variance are calculated 
given a, b, n, and x. 


PREFACE 


The Statistics Library includes this program Manual and 29 prerecorded 
magnetic cards. The programs in this library have been selected because 
of their common use by statisticians and other professionals who need 
rapid solutions to problems involving curve fits, distribution functions 
and general statistics. 


Each program in the Statistics Library is identified by a program name 
and number which also appears on each magnetic card. For example, 
BASIC STATISTICS FOR ONE OR TWO VARIABLES is the name of the 
first program and the number is ST1-01... where ST1 identifies the 
Statistics Library and -01 identifies the first program in the Statistics 
Library. Some programs are in two parts and have two magnetic cards 
which must be run one at atime to obtain different but related 
information. ST1-03-1 identifies the first card of the third program in the 
Statistics Library which is MEANS AND MOMENTS (1). The second card, 
MEANS AND MOMENTS (2), is identified by ST1-03-2. 


The Statistics Library programs are individually described in this manual. 
The general description of a program begins on a left-hand page, 
followed by sample problems, a pictorial of the applicable magnetic 
card, and complete user instructions. The user instructions list all the 
necessary steps to perform the specific program. The first step of the user 
instructions is enter program which should be performed according 

to the procedure on page one of this manual. Each program section also 
includes a program listing and identification of memory registers and 
flags used by the program. The program listing provides a backup 
source to check that the prerecorded program has been correctly read 
by the calculator. 


If you have difficulties in running a program: 
1. Be sure that the user instructions have been followed. 
2. Refer to Maintenance and Service Information in the Owner’s 
Manual to determine the type of trouble you are having. 
3. Write the Consumer Relations Department at: 
Texas Instruments Incorporated 
P.O. Box 22283 
Dallas, Texas 75222 
or call Consumer Relations at 800-527-4980 (toll-free within all 
contiguous United States except Texas) or 800-492-4298 
(toll-free within Texas). If outside contiguous United States, call 
214-238-5461 (We regret that we cannot accept collect calls at 
this number). 


ENTERING A PROGRAM 


Use the following procedure each time you need to store the contents 
of a prerecorded magnetic card in the calculator. Use of this procedure 
is implied by the Enter program instruction in the user instructions of 
each program. 


1. Select the proper magnetic card for the program to be run. 


2. Read side A of the prerecorded 
card as follows: 
iy (Insert card @A ) 
Do not restrict or hold the card 
after it is caught by the drive 
motor. The display is blank until 
the calculator has completed 
reading side A. 


3. Remove the card from the left 
side of the calculator and read 
side B of the prerecorded 
card as follows: 

Fy (Insert Card Be ) 
The disptay is blank until the 
calculator has completed 
reading side B. 


4. Remove card from the left side 
of the calculator and insert into 
upper slot such that the A side 
of the card shows in the window 
above keys A through E. 





5. If the display flashes immediately following step 2 or 3, repeat 
the procedure beginning with step 2. 


CAUTION: Prerecorded magnetic cards may be damaged or altered if 
exposed to dust or foreign materials, permanent magnets, or 
electromagnetic fields (electric motors, power transformers, etc.). 


ST1-01 


BASIC STATISTICS FOR ONE OR TWO VARIABLES 


This program calculates means, standard deviations, and standard 
errors of the mean for one or two variables, and covariance and 
correlation coefficient for two variables. 


The input is a set of data points: 

{(X1,¥1),(X2,y2), oe. (X:,Yi), pee (XnYn) } | 
in which the y; s may be omitted if only one variable is being used. 
The statistics computed are: 


pe 1 n 
x=— x, = mean orx 
MN j=4 
ites Me 
ager 2 yi = mean ofy 
1=1 
n 
2. a He ne 
aad" em n—1 = std.dev.ofxbyn—1method ~ items tr Sam 
n 
p3 } By nx? 
S. = a = std. dev. of x by n method 


= Std. dev. of y by n—1 method 














n—-1 
Dy? — ny? 

os an anna = std. dev. of y byn method 
Sx 

Sx = = std. error of the mean of x byn—1 method weayws 
Vn 

of = Bx Std. error of the mean of x by n method 
Vn 

oS; = Py Std. error of the mean of y by n—1 method 
Vn 

o> = ae Std. error of the mean of y by n method 
Vn 


S X; S y = covariance of x and y 
une by n—1 method 


> Xj és y; | = covariance of x and y by n method 





Pea =a, — correlation coefficient of x and 
ee = ae y 


gr ¥ a 

‘ son m we 

~ BAL AAS ED 
am, ee Wrest ¢ + Aer 
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NOTES: 1. The number of data points must be > 1. 

2. This program may be used with the Linear Regression 
program (ST1-08) to determine additional statistics. Data 
points entered on either program may be used for the other. 
Do not initialize if data from the other program is being used. 


exemple: Calculate x,9,S,, 3,3, 9), Sz 9, 95, S Hie Sx» bry for the 
data in the following table, then change the x and y values fori = 5 
from (23, 17) to (5, 19) and find the new r,,, and S‘,,. 


oO Gon Oo oO 2 @ PhP — 





17 


19 


g>ERES 


PRU ENEEREREREREREREE 


“U 
= 
4) 
¢)) 
ep) 


Display 
i 
iF 
a 
?. 
3. 
3. 

24.38888889 


16.33333333 
10.07403152 
9.497888004 
4.981214711 
4696334268 
3.358010507 
3.165962668 
1.660404904 
1.565444756 
10.85416667 
9.648148148 
2163006968 

8. 

8. 

9. 

9. 
0825567005 
5.006944445 
4.450617284 


Comment 
Initialize 


compute x 

compute y 

compute S, 
compute S‘, 
compute S, 
compute S’, 
compute Sx 
compute St 
compute S- 


compute St 


compute S, , 
compute Sj, 


’ compute r,, 


delete 23 (23, 17) 
delete 17 (23, 17) 
insert 5 (5, 19) 
insert 19 (5, 19) 
compute new r,, 
compute new S,, 
compute new S\, 
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ST1-01 


PROGRAM LISTING 


















000 46 *LBL 0 
001 11 A 46 “a 
002 42 STO 16 *A’ 


003 00 0 
004 01 1 
005 44 SUM 
006 00 0 
007 03 3 
008 40 *x? 
009 44 SUM 
010 00 0 
011 06 6 
012 46 *LBL 
013 34 tan 
014 43 RCL 
015 00 0 
016 00 0 
O17 Sox. 
018 42 STO 
019 00 0 
020 08 8 
021 01 1 
022 44 SUM 
023 00 0 
024 00 0 

25 






























“Denotes 2nd function key 
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PROGRAM LISTING 


46 “LBL 


“Denotes 2nd function key 
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PERMUTATIONS AND COMBINATIONS 


For given n andr where n = r 2 0, the following are calculated. 





Permutations: P(‘) = ery = A-(-31)* «0 h-14+-1) 
Combinations: C(") = Sar 


NOTES: 1. After an n,r pair has been entered, either or both may be 
changed for a new calculation. 
2G) OG) = 
3. Error indications 
a. Forr >n, the display flashes 9.999999999 99. 
b. The display will flash 9.999999999 99 in case of overflow in 
the calculation. 
c. For negative entries, the absolute values will be used. 
d. For non-integer entries, the rounded integer value will be 
used. 


Example: Find the number of permutations of 60 items taken 5 at atime 
and the number of combinations of 120 items taken 5 at a time. 


Enter Press Display Comment 
ea Initialize 
60 [A | 60. n 
9 5. r 
655381440. PCD) 
120 [A | 120. newn 
[ D | 190578024. CCP) 


911-02 












TEXAS INSTRUMENTS 


cO-LLS SN.LVNIGINOOD 8 SN.LVLNINGAd =a» 
@A® PERMUTAT’NS & COMBINAT’NS ST1-02 











Enter Program (A and B) 


linteger(n) | 
linteger(r) | 


a 
C(") 
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PROGRAM LISTING 


084 18 *C’ 
085 42 
086 00 
087 00 
088 43 
089 00 
090 01 
091 42 
092 00 
093 04 
094 75 
095 43 
096 00 
097 02 
098 95 
099 22 
oO 
78 
42 
00 
05 
43 
00 
01 
fa 
07 
00 
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PROGRAM LISTING 





“Denotes 2nd function key 
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MEANS AND MOMENTS 


For a given set of input data, {x1, x2,...x,}, with associated frequencies, 
{fi, fo,...f.,} (for grouped data), the following means, moments, and 
skewness and kurtosis of distribution are calculated. If f, = 1 for alli, the 
calculations are for ungrouped data, otherwise for grouped data. 





Means 
Arithmetic =A =* 2 FX; Harmonic =H = N/ y3 i 
i=1 i=1 i 
! n f. 1/N ; 1 2 1/t 
Geometric’ =G)| = (1 Xi ' Generalized = M(t) = f > f0)') 
=a Yau 
Moments 
mi = an ae 2 ieekan S f(x, — A)4 
Me ity N mo 
m2 = 2 S f\(x; — A)? Kurtosis of distribution = bbe 
N peg as (m2)? 
=~ 3 f(x, —A) Sk f distribution = —™*, 
Ma = ay 2 ewness of distribution = )aa 
N = >» f; n = Number of different groups 


NOTES: 1. The initialization routine assumes a frequency of 1 andt = 1. 

2. The value of t may be changed before any x,’s are entered, if 
changed after any x;, M(t) will be meaningless. 

3. For grouped data if f; = f;_ ,, it need not be entered. 

4. A negative frequency will cause the next x entered to be 
deleted |f| times. 

o. A negative x will cause the geometric and generalized 
means to be invalid. 


Error Indications: 
A frequency ort of 0 will cause the display to flash 9.999999999 99 and 
the old value to be retained. 


Example: The ungrouped (f, = 1) data {x; = 1.1, 2.3, 4.7, 3.5, 1.8, 2.9} 
yields the following: A= 2.716666667, G = 2.451843924, H = 2.187595899, 
m(1) = 2.716666667; m2 = 1.368055556, m3 = .5367592591, 

ma = 3.903019215; kurtosis of distribution = 2.085418497, skewness 

of distribution = .8354470243 


References: Handbook of Mathematical Tables and Formulas, R. S. 
Burington, McGraw-Hill, Fourth Edition, 1965, p. 164. International 
Dictionary of Applied Mathematics, W. F. Frieberger, D. Van Nostrand, 
1960. 


911-03-1, 2 
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PROGRAM LISTING 


000 46 *LBL . 084 00 O 
11 A 065 05-5 

086 55 + 
087 43 RCL 
088 00 0 
089 01 
090 95 
091 
092 
093 
094 
095 
096 
097 
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S1T1-03-1 


PROGRAM LISTING 


a 


46 CEL 
18: °C: a *1/x 
43 RCL 56 “*rtn 
000. 00 
01 00 
55 + 00 
43 00 
00 00 
07 00 
95 = OO 
56 00 
46 00 
19 00 
43 00 
00 00 
08 00 
55 =~ 00 
43 00 
O00 00 
01 00 
95 00 
30 00 
43 00 
00 OO 
03 00 
95 00 
56 00 
46 00 





“Denotes 2nd function key 
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PROGRAM LISTING 


000 46 *LBL 084 09 9 










































































001 11 A 085 55 + 
002 43 RCL 086 43 RCL 
003 01 1 087 O00 0 
004 02 2 088 01 1 
005 40 “x? Gee 75" + 
006 94 +/- 090 03 3 
O07 85 + 091 65 X 
008 53 ( 092 43 RCL 
009 43 RCL O93 01 7 
010 00 O 094 02 2 
011 09 9 095 40 *x? 
O12 55 + 096 40 *x? 
O13 43 ROL 097 95 = 
014 00 O 098 56 “rtn 
015 01 1 099 46 *LBL 
OVW 95 = 100 14 D 
017 56 *rtn 101 11 A 
018 46 *LBL 102 40 *x? 
019 12 B 103.42 STO 


43 104 O1 1 
105 03 3 
106 13 C 
107 46 *LBL 
100 Gr" 7" 


55 - 
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PROGRAM LISTING 


*Denotes 2nd function key 
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RANDOM NUMBER GENERATOR 


This program generates uniformly distributed random numbers 
(0 = U; = 1) by use of the following formula: 


U;., = 7°U,;(Mod 10'°) 
Where: U,, = 1 


Normally distributed random numbers are also generated, using the 
direct method: 


A; rr ( —2in oe COS (27U;. ;) 
Ais, =(—2In U;)'/? sin (27rU; 1) 


Where: U; and U,, , are generated from the above method. 
The final result is: 


Y,=oA,; +m 


which yields normally distributed random numbers with mean m and 
Standard deviation o. The user may input the standard deviation, mean, 
and U,, or the default values are: U; = 0.55, 0 = 1,andm = 0. 


NOTES: 1. Fornormal distribution, the calculator must be in radian mode. 

2. The random numbers return as fixed point 5. 

3. The initial seed is 1.3597878 x 10'°, and is entered when [ A | 
is pressed. 

4. The[e] and Ww are of use only with the printer. [ E | 
initiates continuous printing of U;, and Mw continuous 
printing of Y;. Printing is stopped by pressing and then 
Ea. | 


Example: To generate uniformly distributed random numbers the 
procedure is: | 


Enter Press Display Comment 
LA | 0.00000 Initialize 
0.48405 U; 
0.46815 U2 


Each pressing of [8 } will display a new U;. With a printer, pressing [E | 
once instead of will produce a list of U; which is stopped by 
pressing and then ra. 


To generate normally distributed random numbers the procedure is: 


Enter Press Display Comment 
LA | 0.00000 Initialize 
2 2.00000 mean 
1 ae 1.00000 o 
[ D | 0.91202 ¥a 
[ D | 2.10940 Yo 


ST1-04 


Each pressing of LD _] will display a new Y;. With a printer, pressing 

Wal once instead of [_D ] will produce alist of Y; which is stopped by 
pressing and then _If the mean and o are not entered, the 
default values arem =OQOando = 1. 


References: Theory and Problems of Numerical Analysis, schaum’s 
Outline Series, McGraw-Hill Book Company. Handbook of Mathematical 
Functions, U.S. Department of Commerce, National Bureau of 
Standards, 1964. 
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Printer list of U; 
Printer list of Y; 
Stop step 8 or 9 
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PROGRAM LISTING 


000 46 *LBL 028 O01 1 











































































001 11 A 029 42 STO 
002 01 1 030 00 0 
003 93 « 031 04 4 
004 03 3 032 00 0 
005 05 5 033 42 STO 
006 09 9 034 00 0 
007 07 7 035 05 5 
008 08 8 036 56 *rtn 
009 07 7 037 46 *LBL 
010 08 8 038 12 B 
011 52 EE 039 24 CE 
012 01 1 040 07 7 
013 00 0 041 45 y* 
014 22 INV 042 09 9 
015 52 EE 043 95 = 
016 42 STO 044 65 X 
017 00 0 045 43 RCL 
018 00 0 046 00 0 
019 22 INV 047 00 0 
020 50 ‘stflg | 048 95 = 
021 01 1 049 42 STO 
022 93 « 050 00 0 
023 05 5 051 00 0 
024 05 5 052 55 = 
025 42 STO 053 01 1 
026 00 0 054 


027 06 6 


“Denotes 2nd function key 
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PROGRAM LISTING 




















































































112 00 0 140 42 STO 
113 04 4 141 01 1 
114 93 142 00 0 
lo U5 5 143 30 COS 
116 05 5 144 42 STO 
117 42 STO 145 00 0 
118 00 0 146 07 7 
119 06 6 147 43 RCL 
120 43 RCL 148 00 0 
121 00 0 149 06 6 
122 04 4 150 23 Inx 
123 56 *rtn 151 65 X 
124 46 *LBL 152 02 2 
125 14D 153 94 +/- 
126 60 “if fig 154 95 = 
127 O11 155 30 *Vx 
128 68 *8' 156 42 STO 
129 50 *stflg | 157000 
130 01 1 158 01 1 
131 12 B 159 16 *A’ 
132 59 *7 160 56 *rtn 
133 65 X 161 46 *LBL 
134 43 RCL 162 16 *A’ 
135 00 0 163 65 X 
136 08 8 164 43 RCL 
137 65 X 165 00 0 
138 02 2 166 07 7 





139 167 GO :- 


*Denotes 2nd function key 
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ONE-WAY ANALYSIS OF VARIANCE (ANOVA) 


One way analysis of variance is a technique for testing the differences 
between the population means of k treatment groups, where each group 
i(i = 1, 2,...,k) consists of n; observations x,,(j = 1, 2,...,n\). The 

different groups need not have the same number of observations. 


Formulas used: 


Nj 


Sum of observations in groupi= > x;, 
j=1 


Kk Nj 2 
kK nj (5 2 x: 
Total Sum of Squares =SS= > > x,? —-~———— 





k 
i=1 j=1 | n, 
i=1 
nj k ni 
(2 xu? (2 F xu)? 
Treatment Sum of Squares = TSS = 2 es i 
1= | l ¥ n, 
i=1 
Error Sum of Squares = ESS = SS — TSS 
The value of the F statistic is: 
“a DFe x TSS 
Fr) * BSS 
Where: DF; = Treatment degrees of freedom = k — 1 
k 
DF2 = Error degrees of freedom = & n, — k 
i=1 
Example: NN Xij O73" 
17 Eo 20 tebe IG 18.66666667 #18.55555556 





21/30 45 25 46 19 33 37 33.57142857 85.10204082 
Sh 1S Teh Tae ne 2 12.14285714  43.55102041 


f= 14.11952 
OF: =.2 
Oro 17 
SS = 2692.8 


TSS = 1680.895238 
ESS = 1011.904762 
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We TEXAS INSTRUMENTS 


0-11S JONVIYVA JO TWNV AVM-SNO 2am | | 
/ [Am ONE-WAY ANAL OF VARIANCE ST1-05 
et 








USER INSTRUCTIONS 


mn ee 
eis boon 7 
eS) ell 
P| step 5 mustbe 
|_| performed foreach i) 
yy Display o;? 


8 | Calculate F,DF1,DF2 
iF | Display Treatment DF 
Display Error DF 














ele] | fe] | eel | Bl) BL 
me 
BRRRORRRERORREEE 





911-05 


PROGRAM LISTING 


AY 
56 *rtn 


DISPL 

































































































028 00 0 084 00 0 
029 03 3 085 01 1 
030 40 *x? 086 95 = 
031 94 +/- 087 42 STO 
032 44 SUM 088 01 1 
033 00 0 089 09 9 
034 04 4 090 43 RCL 
og Ut 1 091 00 O 
036 94 +/—- 092 04 4 
037 44 SUM 093 44 SUM 
038 00 0 094 00 0 
039 O1 1 095 06 6 
040 56 *rtn 096 00 0 
041 46 *LBL 097 42 STO 
042 13 C 098 00 0 
043 01 1 099 O01 1 
044 44 SUM 100 42 STO 
045 00 0 101 00 O 
046 00 0 102 03 3 
047 43 RCL 103 42 STO 
048 00 0 104 00 0 
049 01 1 105 04 4 
050 44 SUM 106 43 RCL 
051 00 O TOF O81 
052 02 2 108 08 8 
053 43 RCL 109 SG “rin 
054 00 O 110 46 *LBL 
055 TTT 16:.@ 










ie ne a 
7 S(2x,)2/n, 
ect a Oe Ba I 
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PROGRAM LISTING 


“Denotes 2nd function key 
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t STATISTIC EVALUATION 


(Paired Observations) 


This program evaluates the t statistic 


5 


NH 
t oe 
with n—1 degrees of freedom, to test the difference between the means 


(unknown) of two normally distributed populations, using paired values 
of data from the two populations. 





Given the set of n pairs of x,y,, let 


Ai =Xi— Yi 
—~ 12 
ane Gnd 
£21 (2 a)" 
s.-(5 a (> y) 
et 
Xx; 6 4 7 4 « 
Example: y, 5 2 F 4 2 


DF =4 
enn Ve 
S, = 1.673320053 
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90-LLS G3auIVd) WAS LVLS-} =a i 
/ | @Am t-STAT EVAL (PAIRED ST1-06 
oe es ee a ee ee 






USER INSTRUCTIONS 













Pe mmirembete loa! te ad | ns a 
OE ic dae aay a sek | oe om 
Rea o.| Jonlh oon eis och LAr se otal 
Ss CT CS 7 La ETT 
A CS 
nn | a 
0 LD | a 
compat Si lh a 
Pe | comeawor i cot ti 
Compute ee ee 
[comes eS 
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PROGRAM LISTING 























































































000 56 *rtn 028 95 = 
001 46 *LBL 029 42 STO 
002 11A 030 00 0 
003 42 STO 031 03 3 
004 00 0 032 56 *rtn 
005 01 1 033 46 *LBL 
006 56 *rtn 034.17 °B’ 
007 46 *LBL 035 95 = 
008 12 B 036 94 +/-— 
009 95 = 037 85 + 
010 94 +/- 038 43 RCL 
011 85 + 039 00 0 
012 43 RCL 040 01 1 
013 00 0 041 95 = 
014 01 1 042 94 +/-— 
015 95 = 043 44 SUM 
1016 44 SUM 044 00 0 
017 00 0 1045 02 2 
018 02 2 046 40 *x? 
019 40 *x2 | 047 94 +/- 
020 44 SUM 048 44 SUM 
021 00 0 049 00 0 
022 08 8 050 08 8 
023 43 RCL 051 01 1 
024 00 0 052 94 +/- 
025 03 3 053 44 SUM 
026 4: 054 00 0 

1 055 03 3 





“Denotes 2nd function key 


141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
14 
ive 
153 
154 
1S 
156 
ToT 
158 
199 
160 
161 
162 
163 
164 
165 
166 
167 


PROGRAM LISTING 


140 15 E 

25 GLA 
47 *CMs 
“rset 


86 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
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© 11-07 


t STATISTIC EVALUATION 
(TWO SAMPLE TEST) 


This program evaluates the t statistic 


ay yen 
( | | iy (= ne n. x2 es dV? mit me) 
prety + sniapadna sen TE Nes: TUR comeesats 2S) OR Kmeaiee eee 3) 7 
Ny Ny iy ie ny eh 3 2 
with n, +n, — 2 degrees of freedom, to test the significance of the 
difference (A) between the means (unknown) of two independent random 


samples from two normally distributed populations having the same 
unknown variance. 


Given the sample: 


Xi = [Xi, X2, ..., Xn, 
Fey F 25 se Maid 
then a. + x, 
Nx i=1 
ee Ny 
‘Ai Sd 


Example: x, »74, 82, 105, 110; 124; 111, 126, 128 
¥) 294, 79,80, 9S, 96, 102,115, 117; 129 ,5400 


For A=0, t= .7562582916 
DF = 16. 
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LO-1LS “(A 1d NVS OML) 1WAR LVLS-} =a | 
| @A® t-STAT EVAL (TWO SAMPLE ST1-07 
NE ee ee 


* y A DF INIT 








hee RR eNO 
On a een 
OS a CE 
Ce ae 
eed) 110 00 Sp 
CS ie 
ira: le Seanad 
A ee 
<5 eed 
fos A aed 
‘sot ES i. Meola 
ie Ree 
Pondll Jie to. eae 
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PROGRAM LISTING 


a 084 00 O 
46. LBL. 085 08 8 
11 A 086 
44 087 
OO 088 
01 | 089 
40 090 
44 091 
092 
093 
094 
095 
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PROGRAM LISTING 


112 00 O 
113 

114 

119 

116 

ae 

118 

119 

120 

121 

122 

123 

124 

19 

126 

127 

128 

AS) 

130 

131 

132 

133 

134 

135 00 0 
136 -09 9 
Ml 
138 43 RCL 
139 00 O 





*Denotes 2nd function key 
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1 -08 


—— 


LINEAR REGRESSION 
The linear least-squares fit of input data points (x,y) is calculated using 
the following: 
XXL Sxy 
n 
Slope = m= 
P (2x)? 3x2 
n 
Intercept = b = 2y — Max 


The y’ for entered x and x’ for entered y are calculated as follows: 
y’=mx-+b 
x’ = (y—b)/m 


The coefficient of determination is calculated as: 


r= 234) (or BB 


NOTES: 1. The calculator must be initialized by pressing [NS before 
each new set of data is entered. 

2. The number of data points entered must be greater than 1 
before any of the calculations are attempted. 

3. This program may be used with the Basic Statistics program 
(ST1-01) to determine additional statistics. Data points 
entered on either program may be used for the other. Do not 
initialize if data from the other program is being used. 


Example: 

cin ter ‘21.0 20.1 Zod 
Vi we::!l le Cte ee OG 
r? = .9256678121 

m= 1.22906793 

Db ='—15.40509329 


Fory = 12, x’ = 22.29742931 
For x = 22.3, y’ = 12.00315956 
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intcpt 





Ti] EnerProgam (wand) [ 
25 RE 
Enter data Gs ig 2) 

[| ent siorashicy 
ee Delete data 


oe 
a 
Peep 


a. f? 
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adalat at et else 
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PROGRAM LISTING 


































46 *LBL 
18D’ 
95 = 
94 +/— 
42 STO 
00 0 

02 2 

44 SUM 
00 0 

04 4 

40 *x? 
o4 4+ 
44 SUM 
00 0 

OF 7 

43 RCL 
00 0 
U2-2 

65 X 
43 RCL 
00 0 

OF 

95 = 
44 SUM 
00 O 
OD 

43 RCL 
00 O 


ST1-08 
\ 


PROGRAM LISTING 


O11 

O1 1 

oo “rth 
46 *LBL , 


“Denotes 2nd function key 
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ST1-09 
f 


| ff POWER CURVE FIT 
dive a set of data points: 
1X, Yi), 1= 1, 2, ..., n} 
Where: x; > 0 and y, > O, this program will fit a power curve: 
y = ax” (a > 0) 


The problem may be interpreted as a linear regression problem by using 
the equation in the form: 


Iny=bInx-+tIna 


The computed statistics are: 
1. Regression coefficients 


ae = inv; b> In A) 
n n | 
a (In x;)(In y;) — a ae ny) 
is Re sy, 
(In x;)? ger 


NOTE: nis a positive integer # 1. 
2. Correlation coefficient 


E (In x)(In y,) — & 1 xd In yi)? 


[s(n xy — EHF (in ye — 2 | 


r= 


3. Estimated value x’ for given y 


b 
x'= Vy/a 
4. Estimated value y’ for given x 
‘= ax? 


ii 188 15°48 5) 
So 1.38.03 b6e 106 B4s0 ee 
a = .2901902009 
b = .7759169351 
r2 = 0.921727263 
forx=9 , y’ = 1.596230726 
for y = 2.8, x’ = 18.56904361 














Example: 
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y ‘ 

4 ‘ 
® 
3 & 
2 # 
a 
# 

1 * 
0 x 

0 5 10 15 20 











|| @Am POWER CURVE FIT ST1-09 


a r2 





Xi 


Yi 














ENTER | PRESS DISPLAY 









i | BierProgam hands | 
(1. a eoere 
[teneattorany 





i / e 
x< 


te Calculate correlation 
ions coefficient 


value for 
x’ given y 
y’ given x hatin ee 








“NN 
c 
=) 
= 
<2) 
a 
oP 
cD 
cD 
2) 
np 
= 
jee) 
am 
cD 
co. 


tele | | |e 
Uaeee 
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PROGRAM LISTING 
TLL. 056 46 *LBL 


16 *A’ 
23 
42 
00 
07 
22 
44 
00 
03 
40 
22 
44 
00 
05 
43 


sy -09 


PROGRAM LISTING 


20 *1/x 140 08 8 168 46 *LBL 

42 141 94 169 17 *B’ 

01 142 65 170 55 + 

00 143 43 171 43 RCL 

65 144 00 172 00 0 

53 145 01 173 09 

43 146 85 174 95 

00 147 43 175 35 

04 148 00 176 43 

75 1 149 08 177 00 

43 150 95 = 178 08 

00 151 55 + 179 95 

01 152 43 180 56 

65 153 00 181 46 

43 154 00 182 18 

00 155 95 = 183 45 

03 156 22 184 43 

55 3 157 23 185 00 

A3 158 42 186 08 
159 00 187 65 
160 09 188 43 
161 56 189 00 
162 46 190 09 
163 19 191'\95'= 
164 43 192 56 
165 00 193 46 
166 08 194 15 
167 56 195 43 





“Denotes 2nd function key 
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ST1-10 


| EXPONENTIAL CURVE FIT 


iven a set of data points: 
{(X;, Yi), |= 1, 2,...,n} yi > 0 


This program finds the least-squares fit for an exponential function of 
the form: 


y = ae” (a > 0) 
The problem is computed according to the linear equation 


Iny = Ina + bx 


Statistical outputs: 
1. Coefficients a, b 


Note: nis positive integer + 1. 
2. Correlation coefficient 


| af 
Exiiny as ax2iny 


ai ‘\2 ‘\2 
3x7 | E(iny,) — (lIny;) | 
n a 
3. Estimated value for x’ given y | 


4. Estimated value for y’ given X 


y' — aes 


Example: X, 1.42 1.94 2.63 3.11 se 4.25 





yi .63 9 1.42 Re. 2.43 2.89 


a= .3131935086 
b= .5477927735 
r?= 9810626986 


7524162355 
Fory=2.2, x' =3.558629277 


mn 
O 
oe 
x 
| 
os 
ee 
3 
| 
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STEP 





y’ given x BA TRE Oe 19 
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PROGRAM LISTING 

































001 
002 
003 
004 
005 
006 
007 
008 
009 
010 
011 
012 
013 
014 
015 

| 016 

017 

018 

019 

020 

021 

022 

023 

024 

025 

026 

027 


“Denotes 2nd function key 












10 
43 
00 
00 
96 
46 
11 
47 
25 
96 
46 


000 46 *LBL 


*E! 


028 00 
029 05 
030 43 
031 00 
032 06 
033 44 
034 00 
035 01 
036 49 































038 07 
039 40 
040 44 
041 00 
042 02 
043 43 
044 00 
045 07 
046 44 
047 00 
048 04 
049 01 
050 46 
051 87 
052 44 
053 00 
054 00 
055 10 





037 00 


0 





















96 *rtn 
46 "LBL 
16 *A’ 
23 
42 
OO O 


056 
057 
058 
059 
O60 
061 


















22 
Aa 
00 O 

03 3 

40 *x? 
22 
44 
00 O 
U9 
43 


063 
064 
065 
066 
067 
068 



















070 
O71 
072 






95. = 
20 WX 





REGISTERS 


a ee 
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PROGRAM LISTING 


196 00 0 
197 65 X 
198 53 ( 
199 43 RCL 
200 00 
201 05 
ae 79 
203 43 
204 

205 

206 

207 

208 

209 

210 

211 

212 

213 

214 

215 

216 

217 

218 

2 

220 

221 

222 00 
223 00 





“Denotes 2nd function key 
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ST1-11 


LOGARITHMIC CURVE FIT 
Given a Set of data points: 
ev phe 1, 2,.25R} 
Where: x;>0, this program fits a logarithmic curve: 
y=a-+ binx 


Statistical outputs: 
1. Regression coefficients 


Pay xy; — b2dInx, 


2. Correlation coefficient 


ag 
ay iiInx; — =n ey 
r= 
E(Inx;) ie ariel zy, 24 au 
n n 


3. Estimated value for x’ given y 
x’ = exp b : | 


4. Estimated value for y’ given x 





y’=a-+ binx 
NOTE: nis positive integer + 1. 


Example: 





yi 9.3 10.8 oO Z0.2 36.9 44.4 


a = —13.49953513 
b= 18.91052316 
r2 = 8418448418 


For x = 21, y’ = 44.07397695 
For y = 40, x’ = 16.93001448 
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TEXAS INSTRUMENTS 


LL-LLS Lid SAYND OINHLIYVDOT =a | | 
|_|. @A® LOGARITHMIC CURVE FIT ST1-11 
ae 8 eS eee a 


delete —> x! x= y’ 













USER INSTRUCTIONS 


STEP PROCEDURE ENTER PRESS DISPLAY 







itmctmn © | Pt dT 
mee ee | ue la 
TT eneat tran ts 


> 





set 
‘yeaa 
oe ee 
"Caluaeconeation 
oF eigen og aaa mp 
computeesimsed 
ih ean NR 


fe] | (6) @lel@le| fale 
CORRERORCEREE 
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PROGRAM LISTING 

































084 00 O 
Uso Of 7 
086 22 INV 
087 44 SUM 
088 00 0 
089 04 4 
090 O1 1 
OST 94 +/— 
092 41 GTO 
GS5 Gs “1: 
094 46 *LBL 
095 14 D 
096 43 RCL 
097 00 O 
098 02 2 
oY ko aie 
100 43 RCL 
101 00 O 
2 Or. 1 
103 40 *x? 
ive oot 
105 43 RCL 
106 00 0 
107 00 O 
108 95 = 
TOS gu tex 
110 42 STO 
1h 








ST i+ 





PROGRAM LISTING 
resror [ver [ oaror [wer [oan [ner [ oore [ee 













112 00 0 168 65 X 



























































































































to GO X 169 53 ( 
178 So { 170 43 RCL 
116 43 RCL 177 QD 0 
116 00 O 142.00 
117 04 4 173 75 — 
Le to > 174 43 RCL 
119 43 RCL Ivo uu v 
120 00 O 176 Ua 3 
121 01 1 T77 40 *K* 
1g0.65 A lfe vo. 8 
123 43 RCL 179 43 RCL 
124 00 O 180 00 0 
tes US 3 181 00 O 
tee, oo + 182 54 ) 
127 43,ACL 138320 *“1/x 
| 128 00 O 184 65 X 
129 00 O 185 43 RCL 
130 54 ) 186 O1 1 
131 42 STO 187.074 
ae o) 1 188 40 *x? 
You W711 189 95 = 
134 95 = ‘| 190 42 STO 
135 42 STO 191 O1 1 
136 00 O 102 .-02..2 
137 08 8 193 56° *rtn 
ig0 34 +/— 194 46 *LBL 





139 65 X 195 


*Denotes 2nd function key 
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911-12-1, 2 


MULTIPLE LINEAR REGRESSION 


Using the least-squares method this program (2 cards) derives a linear 
equation of the form: 


Z=ao-7+ aixX + a2y 
for a given set of data points: 
14, Va es P= VS 
According to the normal equations 
2Z; = aon + aidx; + a2dy, 
2X;Zj = AoLX; + arLxX? + a2dxX,y; 
DLYiZ; = aody; + aiLxiy; + a2dy? 
the regression coefficients ao, a1, a2 may be computed. 
Ni — No 
[N&x? — (2x;)?] [nzy? — (Ly,)?] — [n&x.y; — (2x,)(Zy,) ]? 
Where: Ni = [n2&x? — (2x;)?] [nZy;Z; — (Ly;)(=z,) ] 
N2 = [n=xy; — (2x,)(Zy,)] [n=x,z, — (2x,)(2z,) ] 
( 


a2 = 


a1= [NXx;Z; ne (2X; ) (XZ; ) | pl a*|nexy) 4 (x X;) Lyi) | 
oe: n>x? — (=x;,)? 
Beet Bay, Arex) 


ao = 
Nn 


NOTE: nis positive integer + 1 


Card 1 is used when entering data points. If an error is made while 
entering coordinates, the incorrect data may be deleted by entering x; 
and y; in the normal manner. Then, z, should be entered in the display 
and [_D | pressed to delete x, y, and z. When all data points have been 
entered, any sum may be recalled by pressing and the corresponding 
data register number (refer to User Instructions). 


Card 2 must be loaded before calculating the regression coefficients. 
After the coefficients are determined, the estimated value for z,(z’) may 
be computed for a given x; and y;. 


Reference: Computer Methods for Science and Engineering, Robert 
LaFara, Hayden Book Co., 1973 
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1.3 2.6 
1.88 1.24 2.9 
3.1 4.4 3.00 


Example: 












1. Regression line: 
Z = .6462676006 + .4986757882x + 1.242837142y 
2. Regression coefficients: 
ao = .6462676006 


ai = .4986757882 

ao = 1.242637 142 

For: x= 1, y =2.5, z2' = 4.252036243 
Registers: Ro n=4 Fos 3Yi* = 14.692 
Aoi 3%. = 685 Ros SRY? = 14.113 
Ro2 dy; = 7:12 Rar S%.2 = 26.875 
Bee... 22, =. 74.65 Ros Sy;Z,; = 29.899 


Ros py Fe = 13.8525 
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USER INSTRUCTIONS 
Le a oY TS a 
tai ie. 1 ae 
Jie ee 
eens i 
Cw ewe 8 Sk 

[Fate data ne oe a 
2s Oe eo oe 
1 ee Oe sO 
1S... ae CR wk 
ay. be recalled, if desired ee. EQiccit Vas 
es CN OM on a 
i a 
0 ESOT EE xy 
De a es Ce 
Ea SE... ye 
ee ee 
BE scerd-c-hve- sath CMOS | CS EG GS RR 
| eee eee 
i ee 
YE en wsntane cama bc SR MES 
ee er ee 
TT Go cocabiuudhaindi LA Ween EO RG ae ane 
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PROGRAM LISTING 


43 
01 
03 
49 


04 
44 
00 
02 
42 
01 
06 
40 
44 
00 
05 
43 
01 
02 
49 
01 
05 
49 
01 
06 
44 
00 
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PROGRAM LISTING 





*Denotes 2nd function key 
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PROGRAM LISTING 


084 43 RCL 

















































AY 
000 46 *LBL 028 43 RCL 056 20 *1/x 
001 11 A 029 00 0 O57 70 *iferr 
002 43 RCL 030 02 2 058 97 “list 
003 00 0 031 40 *x? O59 | Goi x 
004 00 0 04 060 53 ( 
005 65 X 061 43 RCL 
006 43 RCL 062 01 1 
007 00 0 063 04 4 
008 04 4 064 65 X 
Ue 15 = 065 53 ( 
010 43 RCL 066 43 RCL 
011 00 0 067 00 0 
Ute 017 068 00 0 
013 40 “x? O69 65 X 
014 95 = 070 43 RCL 
015 42 STO 071 00 0 
01 072 08 8 
O73 7a 


43 







085 
086 
087 
088 
089 
090 
091 
092 
093 
094 
095 
096 
















098 
099 
100 
101 
102 
103 

















104 








105 
























106 
107 
108 
109 
110 


O97 


01 


1 


Ob S 


65 
03 
43 
00 
00 
65 
43 
00 
07 
75 
43 
00 
01 
65 
43 
00 
03 
94 
42 


X 


( 
RCL 


S11-12- 2 


PROGRAM LISTING 
































196 G9 9 
197 95 = 
198 41 GTO 
199 14D 
200 00 

201 00 

202 00 

203 00 

204 00 

205 00 

206 00 

207 O00 

208 00 

209 00 

210 00 

211 00 

212 00 

213 00 

214 00 

215 OO 

216 00 

217 00 

218 00 

219 OO 

220 00 

221 00 

222 00 
223 00 








“Denotes 2nd function key 


ad 


911-13-1,2 


QUADRATIC REGRESSION 


The quadratic equation 
Y = a2X* + aiX + ao 
with normal equations 


n 


n n 
> ¥ Santi xe +a & K sheen 
‘ : 


i=] j= p=] 


n n n n 
2 Xi¥i = a2 > x." + a p 3 x + ao 3 X; 
= i= 1 i=1 i= 1 


1 = i= i= 
n n n n 
z Xi7V; = Ge bs “% + al yy xe ao X;? 
= 1 i=1 i=1 >: i=1 


is used to approximate functions of a single variable whose true form Is 
not known, where 


n = number of (x,y) points (n = 3,4,5,...) 
Also, for a given x, this program will calculate 
y(x) = a2x* + aiX + ao 
Where: a2, a1, and ao are determined from the normal equations. 


NOTES: 1. To delete an x,y pair, enter x on{ aA ]andyon . The 
displayed count, n, will be reduced by 1. If you attempt to 
delete an x,y pair before data has been entered, —1 will flash 
inthe display. | 

2. If x; is not entered, it will automatically be assigned the value 
of n when jy; is entered. 


Reference: Computer Methods for Science and Engineering, Robert L. 
LaFara, Hayden Book Co., 1973 


Example: Given the sample: 


x 0 { 2 Oe Eat 6 7 8 
ge, 2 pee eae 7 UO ee) 


find a2, ai, ao, and y(6.2). 


a2 = .2392857143 
ai = —1.505952381 
Go = 3.2 


Thus, y= .2392857143x? — 1.505952381x + 3.2 
and =y(6.2) = 3.061238095 
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Enter Program ST1-13-1 Le GB, ee ea 

2 oe iy ES 
LO Sas 

repeat forall Seccwtes, | Sa aR eee 
2 ne ee 
ee eet 





% 


co 









> 





, 


Display a2 
Calculate y = f(x) 





ST1-13-1 


PROGRAM LISTING 


56 *rtn 

46 *LBL 
15 & 

47 *CMs 
25 CLR 
86 *rset 
46 *LBL 
16 








ST1-13-1 


PROGRAM LISTING 


ae CE 
46 *LBL 
52: Ee 
20. “17x 
O01 1 
94 +/- 
56 “rtn 
00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

OO 

00 

00 

00 

00 

00 

00 





“Denotes 2nd function key 
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$11-13-2 


PROGRAM LISTING 






















































































000 46 *LBL 028 43 RCL 056 20 *1/x 
001 11A 029 00 O 057 70 *iferr 
002 43 RCL 030 02 2 058 97 “list 
003 00 0 031 40 *x? 059 65 X 
004 00 0 032 54 ) 060 53 ( 
005 65 X 033 75 — 061 43 RCL 
006 43 RCL 034 53 ( 062 01 1 
007 00 O 035 43 RCL 
008 02 2 036 00 0 
009 75 — 037 00 O 
010 43 RCL 038 65 X 
011 00 0 039 43 RCL 
012 01 1 040 00 0 
013 40 *x? 041 03 3 
014 95 = 042 75 — 
015 42 STO 043 43 RCL 
1016 01 1 044 00 0 
017 04 4 045 01 1 
018 65 X 046 65 X 
019 53 ( 047 43 RCL 
43 048 00 0 
049 02 2 
050 54 ) 
051 42 STO 
052 01 1 
053 05 5 
054 40 *x? 






055 95 = 







“Denotes 2nd function key 


REGISTERS 










O05 
is oe 






911-13-2 


PROGRAM LISTING 


196 42 STO 






DISPLAY 









































































































TerT°UrTn 
198 02 2 
199 56 *rtn 
200 00 
201 00 
202 00 
203 00 
148 35 + 204 00 
149 43 RCL 205 00 
150 00 O 206 00 
157 OO 0 207 00 
152 95 = 208 00 
183 42.810 209 00 
154 00 O 210 00 
211 00 
56 212 00 
167 46. "LO. 213 00 
158 97 “list 214 00 
159 00 O eho. WO 
160 56 *rtn 216 00 
167 46 "LBL 217 00 
162 12 B 218 00 
163 43 RCL 219 00 
164 01 1 220 00 
165 00 0 221 00 
166 56 *rtn 222 00 





46 *LBL 223 00 


*Denotes 2nd function key 
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ST1-14 


HISTOGRAM CONSTRUCTION 


This program constructs a histogram, given the range (Xmin tO Xmax) and 
the individual observed values (x;) of a variable. Each data value entered 
is checked against the range, and those meeting the criteria 


Xmin = Xj a Xmax 
are placed in the appropriate class interval (cell) by the algorithm: 
Class Interval (Cell No.) = INT (12% —%an_) + 0.5| 


A data value falling outside the specified range is discarded when 
entered and this condition is indicated by a flashing-nines display. The 
maximum allowable data range is 107? to 10°. 


Outputs of the program are: 
n = total number of x, entered within the range specified 
m = mean value of the x; = 2X 


‘<9 
ao? = variance of the x; = t Sx Br oh 


f; = frequency of class interval (cell) j = Mi 
where N; = number of items in cell j 
U; = upper limit of cell j = (Xo Xn | + Xonin 
lf all entries occurin one cell, then the frequency is 1 and the cell number 
which is added to the frequency before display, will appear to be one 


greater than it should be; that is, if all observations fall in cell 7, the 
frequency of cell 7 will be displayed as 8.000. 


EXBINDIOC:? Rein HO, Kear =O 
Wate vaiues: -6,--6,-4,.7-3, —2Z, —3,40.1, 2.3, 4,5 
Number of entries= n= 12 (in the range Xin = X < Xmax) 
Mean = m=—0.5 
Variance = o* = 11.91666667 
Frequency and upper limit versus cell number: 





Cell No. 1 2 3 4 6 
Frequency 

Upper Limit | —5.000 -—4.000 -3.000 -—2.000 -—1.000 0.000 
Cell No. 7 8 9 10 11 12 


Frequency 





UpperLimit;} 1.000 2.000 3.000 4.000 5.000 6.000 
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+\ F 
f 






‘TEXAS INSTRUMENTS 


bl-11S NOILONHYLSNOD WVHDOLSIH =am | 
|| @A® HISTOGRAM CONSTRUCTION ST1-14 
(oon \ meme 2 oe) ee CU 


INIT x (min) x (max) Xi freq 


USER INSTRUCTIONS 








Enter Program (A and B) 


2 Initialize 


Enter lower limit had 
7 | Bisa tl count 
Display variance 

[AT | Desay unpertintot [| ta | | oad 
La), PS a 


Notes: 1) The format of the display in Step 10 is (cell number. frequency). The cell number is incremented by 1 with each press 
of 


2) Step 11 must be performed immediately following step 10 for the cell in question. 


& 


< | x 
3 
7 


EEEROEEEEE 









65 


ST1-14 


PROGRAM LISTING 





“Denotes 2nd function key 


REGISTERS 


09 cell9 14 x(max) 


FLAGS 
Coe ee 
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ST1-14 


PROGRAM LISTING 


140 86 “rset 168 43 Ree . «| 196 01 1 
"Le. . 





*Denotes 2nd function key 
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ST1-15 


NORMAL DISTRIBUTION 


In order to calculate the standard normal distribution, the following 
probability function is used: 


1 





ik} = er} 
0) = on 
Then: Q(x) = Z(x)(bit + beat? + Dal” + Dal* Dsl") + (x) 
wan 1 
Where: t = ng? 4 OX 
e(X)| 4.5 x 10° 
Dp = .2316419 


bi = .319381530 
b2 = —.356563782 
bs = 1.781477937 
b4 = —1.821255978 
bs = 1.330274429 


NOTES: 1.2(+x%) = Z(x), Poche tee Ofx) 
2.1K, 21,0817 
3. Z(x) must be calculated before Q(x) 


Reference: Handbook of Mathematical Functions, National Bureau of 
Standards, 1964. 


Example: For x = 2.02, Z(x) = .0518635767, Q(x) = .0216916245 
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we TEXAS INSTRUMENTS 


GI-LLS NOILNGINLSIG TVNYHON 2am || 
_ | _@A® NORMAL DISTRIBUTION ST1-15 
SS ee ee ee ee 


Q(x) 





USER INSTRUCTIONS 


Enter Program (A and B) 


Enter data 


3 Calculate Q(x) x) 
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ST1-15 


PROGRAM LISTING 


46 *LBL a MOTT 65 43 
11 A 95 = 93 00 
22 INV 42 03 02 
50 *st flg 00 01 45 
01 1 01 09 03 
80 *if pos 56 03 65 
> il 46 08 01 
aa" ty 12 01 93 
50 ‘st flg 25 05 07 
{ 93 03 08 
*LBL 02 76+ 01 
ma? 03 43 04 
01 00 07 
06 02 07 
04 45 09 
01 02 03 
09 65 07 
65 93 75 
43 03 43 
00 5 05 00 
00 06 02 
85 05 45 
01 06 04 
95 = 03 65 
20 07 01 
42 08 93 








“= if 


w 


Oa owawovd-- We x 


w~ 


= A © 
NW ON N FP | ON Cf + K W 


4+NONWOA ODA We XN 





ST1-15 


PROGRAM LISTING 


112 0 

113 vs 
114 05 
119: OS 
116 09 
117 O07 
118 08 
119 65 
120 43 RCL 
121 00 
122 02 
123 45 
124 05 
igo 65 
126 O1 
laf 9S 
128 03 
129 03 
130 00 
131 02 
132 07 
138 04 
134 04 
135 02 
136 09 
137 95 = 

138 65 X 

139 43 RCL 


+ OnNOOaD 


<= NO O 


~*~ 


= < oO 


On FAR FPN DNO OW DW 





“Denotes 2nd function key 
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ST1-16 


CHI-SQUARE DISTRIBUTION 


The chi-square density function is calculated by 


x ‘cme, 
2 a 
pie 


where v = degrees of freedom, 0 < v = 141,x > OandInx = pie 


The following series expansion calculates the cumulative distribution: 


200 In10 
Sh 


P(x|v) = (x) 


Where: ['(v+1) = vI'(v) = v! = v(v—-1)! 
and r(5) eed 


(2) i rk 3 roe eee ramets ay 





NOTES: 1. If the degrees of freedom is entered as a non-integer, it is 

rounded up to an integer for the calculation. 

2. If vy or x is negative or zero, a flashing O is returned. 

3. Av must be entered before an x, and both must be entered 
before P is calculated. 

4. A new x may be entered without changing v. 

o. An x must be entered each time P is to be calculated. 

V 


Example: Forv = 10, r(5) = 24 and if x = 8.5, f(x) = .0969533757, and 
P = .41988 16863 


Reference: An Introduction to Probability and Stochastic Processes; 
James L. Melsa, Andrew P. Sage; Prentice-Hall Electrical Eng. Series, 
1973 
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We TEXAS INSTRUMENTS 


9I-1LS NOILNGILSIG AYVNOS-IHD =am 
@A CHI-SQUARE DISTRIBUTION ST1-16 
CaS: RE, ae 





vp —> E(v/2) 


USER INSTRUCTIONS 


Enter Program (A and B) 
Enter degrees 
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ST1-16 





000 
001 
002 
003 
004 
005 
006 
007 
008 
009 
010 
011 
012 
013 
014 
015 
| 016 
017 
018 
019 
020 
021 
022 
023 
024 
025 
026 
























*Denotes 2nd function key 


A (1... 
ne LT A CES CME 











46 *LBL 
43 RCL 
Ges 

44 SUM 
00 0 

08 8 

43 RCL 
00 0 

02 2 
0° + 
43 RCL 
00 0 

08 8 

65 X 
43 RCL 
00 0 

04 4 

85° + 
42 STO 
00 0 

04 4 

43 RCL 
00 0 

06 6 

95 = 
52 te 
22 INV 
S2 EE 


PROGRAM LISTING 







057 
058 
059 
060 
061 
062 
063 
064 
065 
066 
067 
068 
069 
070 
071 
072 
073 
074 
075 
076 
077 
078 
079 
080 
081 
082 












































056 04 4 
49 *PROD 


00 
03 
46 
42 
43 
00 
04 


0 

3 
"hk 
STO 


io — 


93 
05 
95 
22 
90 
Aq 
59 
30 
49 
00 
03 
43 
00 
03 
96 
46 





085 
086 
087 
088 
089 
090 
091 
092 
093 
094 
095 
096 
097 
098 
099 
100 
101 
102 
103 





























Uoe o2 EE 


o/ 
00 
as 
22 
92 
oF 
09 
22 
80 
10 
90 
10 
42 
00 
07 
95 
01 
42 
00 
03 
02 
95 
42 
00 
01 
42 


“fix 
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PROGRAM LISTING 





“Denotes 2nd function key 


56” 

46 il 

12 B 

22 INV 

52 EE 
INV 
“if pos 


#9 


ort-17 


t DISTRIBUTION 


The t distribution is evaluated by the following series expansion: 


A(t|v) = bs 0 Sin OPOOS P tone? 64 4 ee cos’-2 6 
T 3 OO (E82) 
for v>1 and odd 


A(t|v) = <e for v=1 


A(t|v) = sin 6 42 ens! ee ode g ee -— cos’~? a | 


2 2°4 2-46. ... (v-2) 
for v even 
Where: v is the degrees of freedom 
t >A 
t 
g=tan"' —- 
Vv 
Then: F= Athy +4 aah 


NOTES: 1. Calculations must be done in radian mode. 
Example: For t= 1.341 and v = 15, F = .9000625642 


Reference: Handbook of Mathematical Functions, Abramowitz and 
Stegun, National Bureau of Standards, 1964 
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ST1-17 











x TEXAS INSTRUMENTS 


Zt-t1S NOILNGINLSIG?} sam I | 
[A t DISTRIBUTION ST1-17 
ea a 






DISPLAY 


Slide Deg/Rad switch 
toR 


ag 


ay freedom 
Calculate F 


LF 


ST1-17 





PROGRAM LISTING 
46 *LBL 30 “Vx 





11 A 95 = 





11-17 


PROGRAM LISTING 


112 00 O 140 06 6 
113 03 3 141 55 + 
114: 59 > 142 53 ( 
119. 39 “a 143 43 RCL 
116 95 = 144 00 O 



























117 
118 
119 
120 
121 
hee 
123 
124 
Yeo 
126 
Tes 
128 
V29 
130 
131 
132 
133 
134 
135 
136 
137 
138 
139 
























*Denotes 2nd function key 


10 
96 
46 
93 
01 
42 
00 
06 
51 


*E’ 


145 
146 
147 
148 
149 
150 
Tat 
152 
153 
154 
130 
156 
Ve 
158 
Tae 
160 
161 
162 
163 
164 
165 
166 
167 


























06 6 
oo + 
7 
54 ) 
65 X 
43 RCL 
00 0 

03 3 

33 COS 
40 *x? 
95 = 
42 STO 
00 0 

Go: 5 

44 SUM 
00 0 

04 4 

02 2 

44 SUM 
00 0 

06 6 
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511-18 


F DISTRIBUTION 


The F distribution is most commonly used in terms of the ratio of two 
independent variables with chi-square distributions. 


lf Y and Z are independent random variables distributed by chi-square 
Y/11 _ 
Z/v2 
called the variance ratio. Its distribution is the F distribution with 11 and 
v2 degrees of freedom. 





laws with v1 and v2 degrees of freedom respectively, then F = 


This program calculates Q(x), the probability of an event ‘x >X’’. 
Q(x) = Pr{x>xX} 
Q(x) is the “upper tail’ area on the F distribution curve. 


f(X,v1, 12) 





Q(0) = 1.0 (total area under curve = 1.0) 


This program calculates Q(x) using series expansions, provided that 
either 11 Or v2 is even. 


Conditions: x20, 11 and v2 must be integers > 0 and one must be even. 
Q(x) is calculated according to the formulas: 
1. vi even 


a: Yo {R—K) . v2 (potent —K)? 
Q(x) = K2/1 + 28) 1 te 


Rew? : 


+ eet) ee) 





eS errs 
2. v2 even 
wi 2 
66) = 1-4 0-—Kp i + OR re + ‘ 
uA vi(vit2) ... (v2tvi—4) K © 
O64 ... (p2—2) 
aes 1 
core 5 1 + 11x/v2 
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o1T1-18 


If both 11 and v2 are even, either formula may be used. Choosing the 
smaller of the two will reduce computation time, since fewer terms are 
generated by the series expansion. For example, if v1 = 4 and v2 = 16, 
use the ‘v1 even”’ key Ml to calculate Q(x) more quickly. 


Example: Forv1 =5, v2=8: Q(1.0) = .4748134031, Q(9.5) = .0032423325 
For 11 = 8, v2 = 5:Q(1.0) = 5251865969 











ri TEXAS INSTRUMENTS 


SbL-lLS NOILNGINLSIGS =a =| 
_|_@A« F DISTRIBUTION ST1-18 
See ee eee eae 













STer 


Enter Program (A and B) 
) 


a < 


Inputs (in any order 


NO 


( 
Compute Q(x): 


Se) 


eT ae: 
ee ee 
4 | Ifitis desired to Heo 

recall inputs ian dient 
§_| Tochangeoneormore | 


of the inputs, re-enter 


the desired change 


(step 2) and continue 


to step 3. 





| eleva) fele| [elle 
ete | Et Ey 





ST1-18 
















001 32 
002 43 
003 00 
004 09 
005 85 
006 53 
007 24 
008 75 
009 43 
010 00 
011 06 
012 85 
013 02 
014 54 
015 90 
1016 34 
017 02 
018 54 
019 42 
00 
































46 *LBL 


sin 
RCL 


PROGRAM LISTING 






086 
087 
088 
089 
090 
091 
092 
093 
094 
095 




















































10 Used 
11 


02 x 


084 42 
085 00 0 





04 
94 
85 
01 
04 
45 
43 
01 
00 
04 
65 
93 
01 
85 
93 
93 
41 
32 
46 
34 
54 
00 
95 
06 
46 


STO 





o11-18 


PROGRAM LISTING 





“Denotes 2nd function key 
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ST1-19 


BIVARIATE NORMAL DISTRIBUTION 


WA Ewha ewe roel 
ener, A p* va at") 
Where: uw, = mean of x 
fly = mean of y 
oO, = Standard deviation of x (a, 4 0) 
oy = Standard deviation of y (a, #4 0) 
p =correlation coefficient between x and y (p? < 1) 
Error Indications: +, =0, 0, =0, p? = 1 


Reference: Handbook of Tables for Probability and Statistics, edited by 
W. H. Beyer, The Chemical Rubber Company, 1966 





fixy) re 


Example: uw, = 1,0, = 1.2 
My = 1.3, 0, = .7 
p = 5 
f(2,1) = .0960208596 
f(1,1.3) = .2187813077 
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3 TEXAS INSTRUMENTS 


6L-LLS LSIG IWINHON SLVINVAIG 2am | 
|_| A BIVARIATE NORMAL DIST ST1-19 
ee ee ee ee eee ee es 


Mx ox MY oy p 


USER INSTRUCTIONS 


STEP PROCEDURE ENTER PRESS 





DISPLAY 













ORS OTS eee a a ee eee 
Oe ag ee Na Ce A EE AES 
OSE eee ae et SE ee 
Beg | Pad OT) 
i SR OR 
A A a AER 
a Ak | A RI 
eee) ie eee 
Pe Calculate f(x,y) ei ke f(x,y) 
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ST1-19 


PROGRAM LISTING 


AY 
46 *LBL 43 
11 A 00 













42 03 
00 56 
00 46 
06 1 
46 42 
12 00 
42 04 
00 40 
01 94 
20 85 
43 01 
00 95 
01 22 
56 80 
46 92 
13 42 
42 00 
00 05 
02 20 
43 
00 
04 
96 
46 
16 
42 
*Denotes 2nd function key 
REGISTERS 
A CTR CNS RED CANES: 
eo ee es 





ST1-19 


PROGRAM LISTING 
































1972.43 RCL 
113.00;:0 
114 07 7 
rte 40 *x? 
116 95 
117 55 + 
118 02 2 
119 55 + 
120 43 RCL 
121 00 0 
Tee US 5 
123 95 = 
124 94 +/- 
125 22 INV 
126 23 Inx 
127 55 + 
428 43 ACL 
129 00 0 
130 06 5 
131 30 *Vx 
132 55 + 
133 Cz*2 
134 55 + 
135 Soir 
136 55 + 
137 43 RCL 
138 00 0 
139 


“Denotes 2nd function key 
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911-20 


LOGARITHMIC NORMAL DISTRIBUTION 


ee XV = cr - - ore : 


Where: x is a random variable whose logarithm is normally distributed 
with mean m and variance o?. (x > 0, o? # 0) 


For given m and o”%, the following are computed: 


median = e™ 
mode = e(™~®”) 
mean = emt o/2) 
variance = e™ +?™(e%— 1) 
f(x) 
Reference: Handbook of Probability and Statistics with Tables, 
Burlington and May, McGraw-Hill 

Example: m = 2, 0? = 1.5 


median = 7.389056099 f(.01) = .0000157407 
mode = 1.648721271 f(.5 ) = .058064287 
mean = 15.64263188 f(1.0) = .0858628159 


variance = 851.9412262 









TEXAS INSTRUMENTS 
...dtdddssisitit(wO:..UUtwt(twitd#déC##(ld#dsltlUCee Ut 







1 [enerProvamiaom [|_| +f 
2 | Computemedan——_m —_‘[cal| | 
ae 
C4 | Computemean | —SS~«d TT 
"S| Comput varance «ica | | _—| vance 
i en ee 
am ee 
OW Wi, 


re ira] 
a 
ame RES HH 









STEP 





additional values of x 


S1T1-20 


PROGRAM LISTING 


































































































000 46 *LBL 028 00 0 056 95 = 084 55 + 
001 11 A O29 02.2 057 22 INV 085 43 RCL 
002 42 STO 030 95 = 058 23 Inx 086 00 0 
00 O 031 56 “rtn 059 65 X 087 O01 1 
032 46 *LBL 060 53 ( 088 95 = 
033 13 C 061 43 RCL 089 22 INV 
034 43 RCL 062 00 0 090 23 Inx 
035 00 0 063 03 3 O91 55 + 
036 O01 1 064 75 — 092 43 RCL 
OS? 56 4 065 O1 1 093 00 0 
038 02 2 066 54 ) 094 04 4 
039 95 = 067 95 = Goo oO + 
040 22 INV 068 56 *rtn 096 53 ( 
041 23 Inx 069 46 “LBL 097 02 2 
042 16 *A’ O70. 16. 098 65 X 
043 56 *rtn 071 42 STO O99 59. “sa 
044 46 *LBL 072 00 0 100 65 X 
045 14 D 073 04 4 101 43 RCL 
046 43 RCL 074 23 Inx 
047 00 0 hes il = dam 
048 00 0 076 43 RCL 
049 65 X 077 00 O 
050 02 2 






051 95 = 
Oe: 88° :+ 
053 43 RCL 
054 00 O 
Uon OT 1 
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WEIBULL DISTRIBUTION 


This program can be used to calculate: 
1. f(x1) = mnx® = Vemx) 
Where m>0, n>0O, x>0 
2. Q(x1) =eo™? 
3. X2 fora given Q, Q>0, such that 
em (+ In eins) 
m 


Example: m= .2,n=.9 


f(2) = .1156338734 
Q(2) = .6885184274 
if Q(x2) = .75, X2 = 1.497701748 







TEXAS INSTRUMENTS 










Je npcmmial 
Enter m 


ceo ga we Mert ag 
eR Ce oF RO Pe 
er | 
4 [eompuetey pe ce fim) 
conpueaee) | cet | fom 
| Soe ——— [oy [|] re 







90 




































46 *LBL 
11 A 


42 


STO 


00 O 


06 
ae 
80 
87 
96 
46 
12 
42 
00 
07 
22 
80 
87 
56 
46 
13 
42 
00 
08 
51 
33 
65 


6 
INV 



























PROGRAM LISTING 
















































00 


911-21 


AY 
028 06 084 07 7 
029 65 085 20 *1/x 
030 43 086 95 = 
031 00 087 56 *rtn 
032 07 088 46 *LBL 
033 65 089 87 *1’ 
034 53 090 25 CLR 
035 43 091 20 *1/x 
036 00 092 56 *rtn 
037 08 093 00 
038 45 094 00 
039 53 095 00 
040 43 096 00 
041 00 097 00 
042 07 098 00 
O46 75 = 099 00 
044 01 100 00 
045 95 101 00 
046 56 102 00 
047 46 103 00 
048 33 104 00 
049 43 105 00 
050 00 106 00 
051 08 8 107 00 
052 45 108 00 

43 109 00 

110 00 
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POISSON DISTRIBUTION 


This program calculates the density function f(x) and the cumulative 
distribution P(x). 
sat m*e7 ™ 

x! 


Where:m>0O 
Bee... 


Reference: Handbook of Tables for Probability and Statistics, edited by 
W. H. Beyer, The Chemical Rubber Company, 1966 


Example: m=2.7 
f(0) = .0672055127 


F( 
P( 


8) = .0047075451 
8) = 9980863725 





TEXAS INSTRUMENTS 
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PROGRAM LISTING 


Posmay | xe [owe | ne [osm | xe | om | ner 
00 44 
00 
Pie 
43 
00 
04 
95 
90 
oo 
43 
00 
01 
oe 4 
53 
43 
OO 
00 
85 
01 
54 








65 
43 
00 
03 
95 
42 
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BINOMIAL DISTRIBUTION 


This program evaluates the binomial density function defined by 


f(k:n,p):= (R) p* q®-» oe ee 
o<23 <7 
a=44 72) 
Rt ie, 


With mean: m = np and variance o? = npq 


Where: p = probability of success on each trial 
q = probability of failure 


For agiven k, n, and p, the function f(k;n,p) gives the probability of k 
successes and (n—k) failures, without regard to their order of 
occurrence, inn trials of an experiment. 


Results of the evaluation are: 


f(k;n,p) 
m 
o2 
and the cumulative distribution 
kK 
PUL ee Oe 


Which is the probability of observing k or less successes inn trials, and 
Paya (i) pray” 
k =k’ 


which gives the probability of observing at least k’ successes inn trials. 
For simplicity, f(k; n, p) is written as f(k). 
With k =0 


Fork >0 


Reference: An /ntroduction to Probability Theory and Its Applications, 
William Feller, John Wiley & Sons, 1968 


Example:n=7 p=0.45k=5 


mM. -=g.15 

me? | we Fae 

(5) =.1172214914 
P(5) = .9642938078 


P’(5) = .1529276836 
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TEXAS INSTRUMENTS 


c-lLLS NOILLNGIYLSIG TVINONIG =O | | 






PRESS DISPLAY 





ee es ae 
Pak co 
bl CS 
cs ae 
et 
oo | 
(2nd) | | 
(2nd) | | 
epost 
ae agers 1S 


“oe 

~ > — 
zl S 
te 


Note: All 9’s (flashing) in the display following a data entry, indicates invalid data entry. 
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PROGRAM LISTING 


O1 1 
00 O 
46 *LBL 





PROGRAM LISTING 


911-23 


140 00 O 
141 02 2 


142 
143 
144 
145 
146 
147 
148 
149 
150 
151 
152 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 


94 
85 
01 
oD 
42 
00 
03 
43 
00 
01 
65 
43 
00 
02 


95 = 


42 
00 
04 
06 
46 
14 
43 
00 


04 


65 
43 
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46 *LBL 
80 “if pos 
00 0 
20" 1/x 
56 "rin 
00 

00 

00 

00 

00 

00 

00 

00 

O00 

00 

00 

00 

OO 

00 

00 

00 

OO 

00 

00 

00 

00 

OO 

00 


911-24 


NEGATIVE BINOMIAL DISTRIBUTION 


This program evaluates the negative binomial density function defined by 


honey 

f(k;r, p) = k p'q* Cee eR 
D pe 4 
qh p) 
k =0, 1, 2,... 


With mean: m = “ and variance o? = *. 
Where: p = probability of success at each trial 
q = probability of failure 


For agiven k, r, and p, the function f(k; r,p) gives the probability that 
exactly k failures precede the rth success. 


Results of the evaluation are: 


f(k; r, p) 
m 
oa? 
and the cumulative distribution 
k en 
Pik = 3 ("TE |) pa 
k’=0 


which is the probability of observing k or less successes in r + k 
trials, and 


: et & ood be 
P= = (TK ')erat 


which gives the probability that there are at least k successes in 
r + K trials. 


For simplicity, f(k; r,p) is written as f(k). 
With k = 0 


Fork > 0 
f(ky = VAD Fh a — 4) 


Reference: An /ntroduction to Probability Theory and Its Applications, 
William Feller, John Wiley & Sons, 1968 


Example:r=5 p=09 k=83 

mM = .5555555556 P(3) =0.99497565 
.6172839506 P’(3) = 0.0256915 

f(3) =0.02066715 


| 
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<0 TEXAS INSTRUMENTS 


ve-lLs ISIC 1VINONIG SAILVDSN 2am | 
-] @Am NEGATIVE BINOMIAL DIST ST1-24 
ee Fee 


o? 













Initialize 
Compute variance Pee me 
hie 
a 


aa 
Compute P(k) 


a 
Oh A= aD RS 


Note: All 9’s (flashing) in the display following a data entry indicates invalid data entry. 
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PROGRAM LISTING 


11 Current (k — 1) 
12 Current k 





911-24 


PROGRAM LISTING 


56 *rtn 
46 *LBL 
80 *if pos 
00 0 
20 *1/x 
56: "rin 
OO 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

OO 

00 

00 

00 

OO 

00 

00 

00 

00 

OO 





“Denotes 2nd function key 
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GEOMETRIC DISTRIBUTION 


For a geometric distribution with mean yp, this program will calculate: 
1. p, the probability of an event occurring 





ovis valence o? =e 


p? 
3. 1X be py 
4. the cumulative distribution P(x) = 1—(1—p)**! 
Where:0 =p =. 


Error Indications: A negative mean or x will cause the display to flash 
9.999999999 99 and those values will be ignored. 


Reference: Handbook of Probability and Statistics with Tables, 
Burington and May, McGraw-Hill 


Example:u=4.2 x=83 


p = 1923076923 
go? = 21.84 

f(3) = 1013291726 

P(3) = 5744174749 









TEXAS INSTRUMENTS 


Enter Program (A only 


Compute probability 
Compute P(x P(x) 
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PROGRAM LISTING 


46 *LBL 

11 A 

22 INV 
“if pos 
a he 
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HYPERGEOMETRIC DISTRIBUTION 


This program evaluates the hypergeometric density function for 
given a, b, n: 


a&x)(nCn — x 
pcg , “rr 
Where: a, b, n are positive integers = 69 
a+d = 69 
X Sai n+x = oe 
X 0 OW S,..., GF 
using the relation 
a — a)(p — n) 
feel wie. Fp Ft 
b= 049; 2). oT) 


to find the cumulative distribution 


The mean and variance are found using 
it Lae 
a+D 


s Ay BOR TD 
ie UE St Ry 





Example:a=10b=8n=/7 


f(0) = P(0) = .0002513826 

f(4) = 3695324284 P(4) = .721719457 
f(7) = 0037707391 P(7) = 1. 

mM = 3.888888889 

a? = 1.118373275 
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i TEXAS INSTRUMENTS 


9¢-1LS LSIG DINLSINOSDYSdAH 2am | 
|| @Am HYPERGEOMETRIC DIST ST1-26 
ES a ee ees Gee 








Enter Program (A and B) 
Enter a 














Be — 
Compute mean 


Compute variance 


oes te 


For new a,b go to 2. 


To compute m, o? for 
ee new a,b perform 
| Spee area 


“Error indications after: 


oe 
ae 
we For new n, go to 4. 
ee 





Step 2 ifa<0, a>69, or non-integral 
Step 3 ifb<0, b>69, b non-integral, or a+b>69 
Step 4 ifn<0, n>69, or non-integral 
Step 5 if x<0, x>69, or non-integral 
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PROGRAM LISTING 


46 *LBL no + 084 16 *A’ 
Ve 7’ re O35 42 S10 
43 43 086 00 
00 00 087 05 
00 02 088 O1 
56 10..- 089 94 
46 43 090 42 
18 OO 091 00 
43 04 092 00 
00 54 093 43 
03 ao "x! 094 00 
56 95 095 08 
46 42 096 42 
11 097 O00 
16 098 07 
492 099 46 
00 100 24 
01 101 42 
102 OO 
103 06 
104 
105 
106 
107 
108 
109 
110 
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PROGRAM LISTING 


196 01 1 
TY? 34 ) 
188. 35 = 
18Y 56 “rin 
200 46 *LBL 
autvecl! *B’ 
202 43 RCL 
203 00 
204 01 1 
205 65 X 
206 43 RCL 
207 00 O 
208 04 4 
209 55 + 
270° 18 *C’ 
bit uo = 
fie oc ttn 
213 
214 
215 
216 
air 
218 
219 
220 
221 
222 
223 












a 














© 
© 





io ©) {2 
~ © © © 


©) ©) © © 
— { & © 
ote Me to oa Oo Oo So So 6 


a 
© 





© 
i) 





“Denotes 2nd function key 
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